Phase transition of g-state clock models on heptagonal lattices 
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We study the q-state clock models on heptagonal lattices assigned on a negatively curved sur- 
face. We show that the system exhibits three classes of equilibrium phases; in between ordered 
and disordered phases, an intermediate phase characterized by a diverging susceptibility with no 
magnetic order is observed at every q > 2. The persistence of the third phase for all q is in contrast 
with the disappearance of the counterpart phase in a planar system for small q, which indicates 
the significance of nonvanishing surface- volume ratio that is peculiar in the heptagonal lattice. An- 
alytic arguments based on Ginzburg-Landau theory and generalized Cayley trees make clear that 
the two-stage transition in the present system is attributed to an energy gap of spin-wave excita- 
tions and strong boundary-spin contributions. We further demonstrate that boundary effects breaks 
the mean-field character in the bulk region, which establishes the consistency with results of clock 
models on boundary-free hyperbolic lattices. 

PACS numbers: 75.10.Hk,64.60.Cn,02.40.Ky 



I. INTRODUCTION 



The role of geometry has continued drawing attention 
in statistical physics. A curved surface, for example, has 
been a useful test ground to study ergodicity [l|, and 
curved nanoscale carbon structures have been expected 
to possess interesting elastic and magnetic properties [j| ■ 
Recently, rapid development of soft material sciences also 
requires a precise understanding of physics on a curved 
surface in terms of geometric interactions One im- 
mediate question from the statistical-physical viewpoint 
is how phase transitions occur on such a curved sur- 
face since they in general depend on geometrical factors. 
In particular, a negative Gaussian curvature yielding a 
saddle-like hyperbolic surface has been more commonly 
studied in critical phenomena than a positive one since 
a positive curvature tends to make a closed surface so 
that it is hard to extend the system size while keeping 
the magnitude of the curvature constant. In a nega- 
tively curved surface, the length scale grows only log- 
arithmically with the surface area, and thus one could 
expect a mean-ficld-like critical behavior in many sys- 
tems. Whereas this expectation was proven true for the 
bulk of the Ising spin system 0, Q , the XY spin model 
has no local order at finite temperatures [|| . This lack of 
order in the XY model is attributed to the gaplcss spin- 
wave excitations that can arise from the boundary at any 
finite temperature T . This argument is based on the fact 
that a negatively curved surface contains a huge amount 
of boundary points: that is, for a negatively curved sur- 
face, the ratio of surface area to perimeter (which is the 
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two-dimensional example of the so-called surface- volume 
ratio in general dimension) remain nonvanishing even in 
the large-system limit. Since it was pointed out that a 
system may have a novel behavior due to the presence 
of a nonvanishing boundary |7j, there have been ongoing 
studies to clarify this issue (ill, H, i, [H [ll[ . While the 
boundary effects can be sometimes excluded, for exam- 
ple, by using a periodic boundary condition [l2| or by 
mathematical abstractions [!, [H, [id], HH , it is often cru- 
cial to understand how a boundary affects the physical 
properties since it may give the most important contri- 
bution to an observed behavior as will be explained in 
this work. 

The complete difference between the Ising and the XY 
models with respect to the presence or absence of the or- 
dered phase motivates us to study the g-state clock model 
on a negatively curved surface. The g-state clock model 
is equivalent to the Ising model for q — 2 and approaches 
to the XY model for q — ► oo. Thereby one can ob- 
tain a better understanding on how the phase structure 
changes in between with varying q. In this paper, we 
present the following findings: first, the critical tempera- 
ture T c is indeed proportional to the energy gap to excite 
the spin fluctuations, second, we report an intermedi- 
ate phase with a diverging susceptibility between the or- 
dered and disordered phases. While it corresponds to the 
quasiliquid phase in the planar case, an interesting differ- 
ence is that this intermediate region characterized by the 
vanishing order parameter and diverging susceptibility is 
observable at every q > 2 on the curved structure. This 
point will be further discussed by studying the Cayley 
tree analytically. 

This work is organized as follows: in Sec. [Ill we explain 
the construction of our lattice for describing a negatively 
curved surface, and introduce the g-state clock model on 
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FIG. 1: Schematic view of a heptagonal lattice with a level 
I — 3, projected on the Poincare disk. 

top of it. The results will be presented and discussed in 
Sec. Mil We then summarize this work in Sec. |lVl 

II. CLOCK MODEL IN HYPERBOLIC LATTICE 

A Schlafii symbol {k,w} means a tessellation that w 
regular fc-sided polygons meet at each vertex. Satisfy- 
ing (k — 2)(w — 2) > 4, every pair of {k,w} results in 
a negatively curved surface, yielding a hyperbolic tes- 
sellation [Ty] . Each hyperbolic tessellation gives a re- 
sulting lattice structure, which will be generally called 
a hyperbolic lattice. In this work, we construct one 
type of hyperbolic lattices, i.e., a heptagonal lattice de- 
noted as {k, w} = {7, 3}, in a concentric way as depicted 
in Fig. Q] We start with the zeroth layer, a point in 
the middle of the Poincare disk and surround it 

by three heptagons. Then the newly added 15 points 
constitute the first layer. Likewise, attaching 12 hep- 
tagons all the way around the first layer adds 45 more 
points, which make the second layer, and so on. A 
heptagonal lattice of a level I means that it is made 
up to the lih layer, and its system size is then given 
by N(l) = 1 + f 5 J2 l j=1 [(^y - (3^I)i]. As N(l) 
increases exponentially with I, the surface- volume ratio 
docs not vanish even in the large-size limit. 

An important consequence of the non-vanishing 
surface- volume ratio is an enhancement of boundary ef- 
fects that exceeds the bulk-spin contributions. Some- 
times only the bulk properties are studied by restricting 
ourselves to a distance less than xl from the zeroth layer 
with a constant < x < 1. However, one should remem- 
ber that the system would not be properly described by 
the bulk part since its fraction eventually vanishes: sup- 
pose that N(l) ~ e zl for some curvature-dependent con- 
stant z. The bulk fraction is then e xzl /e zl = e~( 1-a: ) 2Z , 
which exponentially decreases as I grows. This is why 
the boundary-spin contribution plays a dominant role in 
determining the physical properties of the whole system. 

By the g-state clock model, we mean a spin system 
described by the following Hamiltonian: 

H = -J^Sj ■ sj h- 3j 

(ij) 1 



= -J^2cos(9 t - 6j) - ^/icos(9 4 , (1) 

(ij) i 

where each spin Sj can have one of q possible angles, 
9i = 2irrii/q with rii = 0, 1, . . . , q — 1, and h is a mag- 
netic field along the direction for 9 = with a magnitude 
h. The summation is over the nearest neighbors, and the 
coupling constant J > is the strength of the ferromag- 
netic interaction. As mentioned above, q = 2 and q = oo 
correspond to the Ising and XY models, respectively. In 
addition, the case of q = 3 is equivalent to the three- 
state Potts model The case of q = 4 has the same 
universality class as the Ising system since the partition 
function of the four-state clock model at temperature T 
is formally isomorphic to that of two uncoupled Ising sys- 
tems at T/2 

In the planar case, the g-state clock model for h = 
generally has three phases in the q~T plane (2^1 ■ Two 
among the three are ordered and disordered phases as in 
the Ising model. From the existence of the Kosterlitz- 
Thouless (KT) phase in the XY limit [U , one can argue 
that the third, quasiliquid phase e merg es for q > 4 in 
the intermediate temperature range [22j. The low tran- 
sition point where the ordered phase vanishes is roughly 
described by T\ oc 1/q 2 , as explained in the Villain ap- 
proximation [22|, . On the other hand, the high tran- 
sition point, where disordered phase begins, remains al- 
most constant around T2 = Tkt — 0.89J/fcs for q > 8, 
where fee is the Boltzmann constant [2(| . 

With a constant negative curvature, as shown in the 
next section, some of these behaviors still look qualita- 
tively similar. Specifically, the lower transition point is 
roughly proportional to 1/q 2 whereas the higher one does 
not change much as q increases. However, there also ex- 
ist clear differences in that the intermediate phase be- 
tween these two temperatures is created by a very dif- 
ferent mechanism discussed later, and is present at every 
q>2. 



III. RESULTS 

A. Ginzburg-Landau theory for homogeneous 
lattice without boundary 

Phase transitions on a curved surface can be very dif- 
ferent whether a boundary of the system is considered 
or not. As our numerical experiments include both of 
the curvature and boundary effects, we will first consider 
only the curvature effects in this part, in order to high- 
light the boundary effects more clearly. 

Suppose the g-state clock model is in a continuum 
limit. Phcnomenologically one may write a dimensionless 
free energy F of this system in the ordered phase [24j as 
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where p is a position vector (r) rescaled by a spe- 
cific length scale, £, so that |p| = |r|/£. In Eq. ((2]), 
t) = t)| exp [i(f>{p, t)] is a complex order param- 
eter, h(p, t) is a dimensionless magnetic field represented 
as a complex number, and v is a positive constant. Func- 
tional differentiation of Eq. with respect to ip* yields 
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Assuming the free-energy minimum, SF/Sip* = 0, we 
differentiate Eq. Q with respect to h to find an equa- 
tion for the two-point correlation function, G(p,p') = 
d^(p)/dh(p'): 
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W~ 2 G(p,p') = S(p,p'). (4) 



For a translationally invariant system, we may set p' = 
without loss of generality. Let us take a sufficiently small 
v so that this system has ground states with \ip\ > [25| . 
Then one finds « 1 for small v from which it follows 



(V 2 -l)G(p) = -<5(p), 



(5) 



where p = \p\. 

We now impose negative Gaussian curvature to the 
underlying surface of the model. On a hyperbolic surface, 
the Laplacian operator is replaced by A written as [2(| 



A = 
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where the approximation can be taken due to the expo- 
nential increase of sinh p. Then Eq. (J5]) is reduced to 



dp 2 



+ ^-1 



(6) 



This equation is solved yielding G(p) ~ e -p( 1 +V5)/2 ^ 
Note that the correlation function basically behaves like a 
one-dimensional case as G(r) ~ e~ r ^ Hzl, where £ serves 
as the correlation length of the system. Such an exponen- 
tial decay in G(p) on the hyperbolic surface apparently 
suggests the absence of order-disorder phase transition 
in the current system. This is, however, not the case. 
A noteworthy point is that the number of spins N s (r) 
within a distance r also increases exponentially in a hy- 
perbolic surface, which results in divergence of the mag- 
netic susceptibility defined by 
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FIG. 2: (Color online) Phase diagram of the q-state clock 
model on the heptagonal lattice as shown in Fig. [T] We de- 
fine two transition temperatures T c and T„ so that there exist 
the ordered phase below T c and the normal disordered phase 
above T 3 . The intermediated phase is characterized by a di- 
verging susceptibility with no magnetic order. The dotted 
lines mean extrapolated behaviors of the transition tempera- 
tures to the XY- model limit (q — oo). 



where the summation is over every possible pair 
of spins. Even if G(r) shows an exponential decay, sus- 
ceptibility x is arj lc to diverge at finite T by satisfying 
N s (r)G(r) > 1, from which the critical temperature T c 
can be located Hf| . 

In addition to the curvature effects mentioned above, 
we should take note of the effects of strong boundary- 
spin contributions that are inherent to the present sys- 
tem. Notice that in passing from Eq. j2]) to Eq. we 
have discarded a surface term. As mentioned already, 
however, the boundary effect cannot be neglected in any 
physically realizable system with a constant negative cur- 
vature. Hence, the present system involving the bound- 
ary effects will exhibit distinct pro perties from the mean- 
field character observed in Ref . [14| wherein the boundary 
effects are artificially excluded. We also note that our dis- 
cussion in the previous paragraph supports the validity 
of the mean- field description in the boundary- free system 
since the correlation function decays so fast at T c [271 ]. 



B. Estimation of the Lower Transition 
Temperature T c 

In Fig. [2l we propose a phase diagram of the clock 
model on the physically realizable hyperbolic lattice in- 
troduced in Sec. [TXJ In this diagram, we define T c as the 
temperature above which the magnetic order vanishes. 
Apart from the ordered and disordered phases, we can 
identify the third intermediate one, which is also disor- 
dered but exhibits a diverging susceptibility. Therefore, 
we specify one more transition temperature denoted as 
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T s , above which the susceptibility divergence disappears 
and the normal disordered phase begins. 

In order to obtain the phase diagram, we employ the 
parallel tempering method [29| and measure the magnetic 
order parameter 



<H> = 



where (■ ■ •) represents the thermal average. 
Binder's fourth-order cumulant [301, 



From 
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for different I, we can locate a unique crossing point for 
each q [Figs. [3Ja)-[3]^c)]. This determines the lower tran- 
sition temperature T c as a function of q. 

Figure [3jd) shows the dependence of T c on q. T c is 
found to rapidly decrease to zero as q grows larger. A 
striking observation is that T c is determined by the typi- 
cal energy scale AE to rotate a spin in the fully ordered 
ground state. AE is roughly given by 



AE oc 1 — cos 



= sin 



(8) 



in units of J/ks [sec Eq. (JTJ], being proportional to T c 
for each q as clearly shown in Fig. [3jd). In addition, 
Eq. {8} can be approximated by T c oc 1/q 2 for large q, 
which is analogous to the planar case. More interesting 
is the fact that the relation of T c oc AE captures the 
exact relation, T c (q = 4) = \T c (q = 2), mentioned in the 
previous section. These results are consistent with the 
interpretation that the spin- wave excitation breaks every 
magnetic order in the XY model @; in fact, Eq. ((5J) leads 
to AE = in the limit of q — > oo, and thus T c = 0. 

It is worthy to mention the significant contribution of 
boundary spins to the determination of T c ; this is caused 
by the fact that the actual magnitude of AE depends 
on the number of neighbors. Since boundary spins have 
fewer neighbors, the proportionality constant in Eq. ([5J 
takes a smaller value than those of bulk spins so that 
their orientation will be strongly disturbed by thermal 
fluctuations. We also comment that the spin-wave ex- 
citation is observable in the hyperbolic lattice without 
boundary since it is the basic excitation mode. In the 
latter system, however, the excitation is not sufficient to 
destroy the ordered phase but arises as a separate peak 
in specific heat at q > 4 [Til ]. 



C. Two distinct scaling relations around T c and T s 

We next evaluate critical exponents of the transition 
by employing the finite-size scaling analysis. As pointed 
out in Ref. [7|, distribution functions of \m\ for the hep- 
tagonal lattice deviate from the Gaussian distributions 
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FIG. 3: (Color online) Binder's fourth-order cumulant for the 
g-state clock model on heptagonal lattices with (a) q — 3, (b) 
q — 6, and (c) q = 12. (d) Transition temperatures between 
the ordered and disordered phases, T c , compared to sin 2 (71-/(7). 



(Fig. 2]). Since the idea of the fourth-order cumulant 
assumes a Gaussian peak shape [30j], a direct scaling of 
the cumulant will give a different value from the actual 
correlation-length exponent estimated from the order pa- 
rameter In order to find an appropriate estimate, 
therefore, we perform at each q the scaling analysis for 
(\m\) based on the scaling hypothesis: 



(|m|) oc N~ / 9 f (\T - T^N 1 ^ 



(9) 



In the present case, we choose N instead of / as a proper 
scaling variable, as N gives much better scaling collapse 
at T c than /. Although the finite-size scaling of the 
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FIG. 4: Distributions of the magnetic order parameter for 
Ising spin systems (q = 2) (a) in a heptagonal lattice and 
(b) in a plane. To compare these two cases, we make both 
systems have sizes of N ~ O(10 3 ), and set T « 0.8 T c to 
observe low-temperature regions. A narrow Gaussian peak is 
clearly shown in the planar case while a longer tail is observed 
at low \m\ in the heptagonal case 0]. 



Binder's cumulant with v fails due to the non-Gaussian 
nature of the magnetization distribution, T c 's estimated 
from the crossing of Un and from Eq. ([9|) are almost 
identical. Figure [5] shows the resulting scaling plots and 
estimated critical exponents as functions of q. While v 
appears to be relatively constant at q > 6, (3 tends to 
decrease to zero, suggesting that that every q-state clock 
model belongs to a different universality class, apart from 
the exact equivalence between q = 2 and 4. 

Measuring the magnetic susceptibility x = ^((\ m \ 2 ) ~ 
( m |) )/k B T usually gives another way to estimate v with 
a similar scaling hypothesis, 



(10) 



This yields consistent results with the above ones for q < 
4, and confirms the results in Ref. [i| for q = 2 and 4. 
However, we find Eq. (JTUJ) inapplicable at q > 4 to obtain 
critical indices since the susceptibility begins to diverge 
at a temperature T s , much higher than T c . In contrast, 
the length scale I successfully works as a scaling variable 
(Fig. [5]). Henceforth, we should employ the following 
alternative scaling hypothesis around T s , 



X <xl^'^g(\T-T s \ I 1 '"), 



(11) 



which locates the phase-separation point T s as depicted 
in Fig. [21 In a usual e?-dimensional lattice, there exists a 
trivial relationship between exponents found in Eqs. (|10p 
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FIG. 5: (Color online) Finite-size scaling using Eq. (|9]) for 
(a) q = 4 and (b) q = 16. Here the parameters are chosen 
as p/v = 0.167, 1/v = 0.285, and T c = 0.63 for q = 4, and 
P/v = 0.0075, 1/v = 0.275, and T c = 0.05 for q = 16. By this 
way, we estimate behaviors of critical indices (c) P and (d) v 
as q varies. 



and (fTTj) . derived from N ~ l d . In absence of such a re- 
lation between N and I, it is rather nontrivial to observe 
these different scalings in a single system at different tem- 
peratures. A similar change in the scaling variable across 
two transitions is also found in percolation phenomena on 
hyperbolic lattices [Toj . 

It is noticeable that a diverging susceptibility at finite 
T s appears to be the counterpart of the susceptibility 
divergence at Tkt in the planar XY model. More inter- 
estingly, the higher transition temperature T s , separate 
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(Color online) Susceptibility scaling with / in case of 
(a) The crossing point at T s = 0.6 with ■y 1 /v = 2, and 



(b) scaling collapse with 1/v = 1.5. 



from T c , exists for all q for the heptagonal lattice, whereas 
the quasiliquid phase in the planar case does not appear 
with q < 4. To look into its origin, we below examine the 
clock model on the Cayley tree. 



D. Comparison with Cayley tree 

The Cayley tree is a special type of hyperbolic lattices 
containing no loops, which often allows exact calculations 
as a useful guidance. In order to understand the existence 
of the intermediate phase, we extend the results for the 
Ising model ( q = 2) on the Cayley tree, presented in 
fiefs. HH and |32|. to general q. 

Let us consider a branching number of B = 2, i.e., 
a binary tree with n generations, where a root node is 
denoted as the zeroth generation. The total number of 
nodes are N n = 2 n+1 — 1. Let Z„ denote the partition 



function of this branch, restricted that the root node has 
a phase variable as 9. The complete partition function 

would be then Z n = 

2zr 



AS) 

1 



J2e Zn 1 , where the summation 



runs over 6 = 0,^-,^-,..., "' w g 11 . A tree with (n + 1) 
generations can be generated by attaching two trees with 
n generations to a single node, which has a certain angle 
9. Since two n-generation trees are totally independent, 
one can write the following recursion relation, 



7W _ /3?icos6> 



y"^(0') e /"cos(0- 



(12) 



where (3 = (^T)" 1 and the magnetic field h is assumed 
to be in parallel with 9 = 0. 

By differentiating Eq. (fl"2)) and taking the limit as h — > 
(see Appendix we find the magnetization of the n- 
generation tree with broken symmetry as follows: 
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where R = 
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Following 



the argument in Ref. [33j , we remark that the correlation 
between a pair of spins, separated by the distance r, is 
given as R r . Since R < 1 in general, the magnetization 
(m) n goes to zero as n — > 00. One may also consider 
the free-energy cost of forming a spin cluster on a sub- 
branch of this tree. Since a single bond divides the whole 
tree into two regions, the energy cost at the interface is 
5E = 2 [1 — cos(27r/q)], basically constant regardless of 
the cluster size. At any finite temperatures, the entropy 
gain 5S, by forming a spin cluster, will thus dominate 
the free-energy change, readily breaking the magnetic or- 
der (see Ref. [Hf for a typical spin configuration). Yet 
one should note that this large-system limit can be quite 
subtle [H. 

The second-order derivative of Eq. (fT2")) leads to (see 
Appendix [C| 
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where S = ^ 6 



This formula recovers the Ising case with q = 2, where Xn diverges at (3 S J 
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J/k B T s = ln(l + a/2) w 0.8814 [11]. Also in gen- 
eral, Eq. |H]) diverges at i? = l/\/2. If g -> 
oo, we may rewrite the summations in R as inte- 



grals so that R 
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Ii(j3J)/Io(/3J), where I n (x) is the modified Besscl func- 
tion of the first kind. A numerical solution then gives 
(3 S J as 2.0582. Since Xn{R) is basically the same at 
any q and R is always a monotonic function of tem- 
perature, the divergence should be also qualitatively the 
same as in g = 2. That is, susceptibility diverges as 
Xn/n ~ Qo + &i(T — T s )n with some constants ao and 

ai im. 

The susceptibility divergence can be explained by the 
presence of boundaries [33J. From the viewpoint of 
boundary spins, which dominate the overall property, 
the effective number of generations appears as n' ~ 2n 
since it is the maximum possible distance in this tree. 
Therefore, the effective branching number for a bound- 
ary spin amounts to B ~ \/2 so that B n = N ~ 2" (see 
Ref. 3] for a general discussion). According to Eq. (0, 
the contribution of each boundary spin to susceptibil- 
ity is roughly C = & ffl . Since the number of 
boundary spins is proportional to the system size N, we 
find the lower bound of susceptibility that \ n > f3C, 
and expectedly this will make the most dominant term. 
At R = l/y/2, we have BR ~ 1, which means that 
C ~ n 0(1). Note that the summation is limited by 
the number of generations, n. In other words, the sus- 
ceptibility diverges with \n oc n at R = 1/ v2. 

Recalling differences between with and without loops 
in percolation phenomena [ic| |. we may expect only a 
qualitative understanding for the heptagonal lattice from 
studying the Cayley tree rather than a quantitative 
agreement. Although the presence of closed loops will 
presumably alter the results described above, the essen- 
tial parts of these arguments could be conveyed to our 
heptagonal lattice. That is, the susceptibility divergence 
at T s should be attributed to the exponential growth of 
N(l). This is markedly different from the case in regu- 
lar lattices, where the susceptibility divergence is due to 
divergence in the correlation length. In particular, the 
correlations among boundary spins play the most impor- 
tant role at this point. Nonetheless, the correlation func- 
tion does not have to decay algebraically yet, which is a 
possible reason that Binder's cumulant does not detect 
T s . One cannot observe the algebraic decay until reach- 
ing T c (< T s ). Around that point, the hyperbolic lattice 



begins to manifest itself more as a surface. In contrary 
to the tree case above, for example, the energy cost at 
a domain wall increases roughly logarithmically with the 
cluster size (Tol | . opening the possibility for T c to be fi- 
nite. As a consequence, we observe these three phases in 
general: an ordered phase, a disordered phase but having 
a diverging susceptibility, and a normal disordered phase 
with a finite susceptibility 

IV. SUMMARY 

We investigated the g-state clock model on the heptag- 
onal lattice, and found that the spin-wave excitation is 
relevant in the order-disorder transition in this system. 
In the planar g-state clock model, one could expect one 
additional quasiliquid phase, and thus two phase transi- 
tions for q > 4. The lower transition defines the line be- 
tween true- and quasi-long-range order, and the higher 
one defines where the quasi- long-range order vanishes. 
If we only introduce the curvature effect but without 
the finite surface-volume ratio, the quasi-long-range or- 
der becomes a genuine order and the higher transition is 
of the mean-field type since fluctuation decays exponen- 
tially (see Sec IIII A|) . However, the presence of a bound- 
ary cannot be neglected, which breaks the mean-field pic- 
ture, and the spin-wave excitation appears to be crucial 
in establishing the ordered phase. In the limit of q — > oo, 
the excitation becomes gapless so that the transition tem- 
perature approaches zero. In addition, the susceptibility 
begins to diverge at a higher temperature, indicating a 
similar phenomenon to the KT transition with a diverg- 
ing susceptibility. By analyzing the clock model on the 
Cayley tree, we suggest that the hyperbolic nature of the 
underlying lattice structure makes the third phase ob- 
servable for every q > 2. 
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APPENDIX A: MAGNETIZATION IN CAYLEY TREE 

If we take the limit of h — > 0, Eq. (fT2"|) leads to 



z (e) _ Z {9y 
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3 /3Jcos(0-0') 
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since Zn^ = Zn °' ) by symmetry. It is straightforward to see that 



1 2JV„ 
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since 



(A2) 
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Note that Eq. (|A2|) is an analytic function at any T [3j|. As to derivatives, one finds the following equations by 
differentiating Eq. (fT2|) with respect to j3h: 
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where 



We then take the zero- field limit, h — > 0. By mathematical induction (see Appendix |B|) . one can show 
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Henceforth, by Eqs. (|A1[) and (|A5|) . we can rewrite Eq. (|A4[) for a restricted ensemble with 9 = as follows, 
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which directly leads to Eq. ([13 
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APPENDIX B: MATHEMATICAL INDUCTION 



Let us assume that Eq. (|A5|) holds true, as it does for n = 0, 



<9(/%) 



h=0 



h=0 



COS 9. 



Then for general n, this assumption yields the following relation: 
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Here we note the following identity: 
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where 9" = 9 — 9' and the last equality is due to the fact that sin 9" is an odd function. Therefore, we substitute 
Eqs. (TM|) and (TB2]) into Eq. (TBI]) and then obtain 
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which confirms Eq. (|A5[) for any n > 



cosy, 



APPENDIX C: SUSCEPTIBILITY IN CAYLEY TREE 



For describing susceptibility, we again differentiate Eq. (|A4[) to get 



d{{3h) 2 d(ph) 
In the zero-field limit, we have the following: 
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To simplify the last term, we sum up both sides over 9 and find 



1 ag 



+ 2- 



1 d 2 zV>') 0j cos(e _ 6 >) 



0' ^W) d((3h) 2 



T,e> e 



/3 J cos 8' 



E 



i ^ 

^ dm 2 



S{{\+ AR) 



(o) 



i az, 



2i? 2 



1 0Z 



(0) 



^E 



1 d 2 zk e) 



(C2) 



where 5 = cos 2 0. Now Eq. (|C2|) describes the full ensemble without breaking symmetry, which is valid above 
criticality. The terms inside the curly brackets can be explicitly written by using Eq. (|13[) . Solving this recursion 
relation with the first term as 

we obtain Eq. (TTJ} as the susceptibility for the n-generation tree. 
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